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Detailed  observations  of  a  naturally  occurring  shear 
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Abstract.  Simultaneous  profiles  of  microstructure,  horizontal  velocity,  and  acoustic 
backscatter  allow  one  of  the  most  complete  descriptions  of  a  naturaUy  occurring  shear 
instability  to  date.  Shear  increased  rs^iidly  after  passing  through  a  lateral  constricticm 
which  formed  a  hydraulic  control.  A  kilometer-long  set  of  20-m-tall  billows  grew  on  a 
middepth  density  interface  where  the  Richardson  number  fell  below  0.25.  The  velocity 
interface  thickened  steadily  after  the  billows  formed,  consistent  with  rapid  momentum 
mixing  across  a  shear  layer  with  a  Reynolds  number  of  3  x  10^.  The  billows  generated  large 
density  overturns  and  dissipation  rates  greater  than  10~’  W  kg~^,  even  within  the  first  large 
overturn,  indicating  that  these  structures  were  fully  turbulent  early  in  their  development. 
As  the  billows  grew,  a  well-mixed  layer  developed  at  the  interface  and  survived  as  an 
actively  turbulent  layer  for  up  to  6  buoyancy  periods,  3  times  longer  than  in  laboratory 
studies  at  low  Reynolds  number.  Variations  in  the  mean  density  of  the  billows  lead  us  to 
infer  that  the  vertical  offset  of  the  velocity  and  density  interfaces  vuied  with  time  where 
the  billows  first  formed.  With  data  from  the  large  overturns  within  the  shear  layer,  we 
find  e/vN^  «  3  x  10'*,  an  average  root-mean-square  overturn  scale  (Lrnu)  of  2.6  m,  and  a 
buoyancy  scale  (L^)  of  2.7  m.  Despite  having  sampled  the  biUows  at  varying  stages  of  their 
evolution,  we  find  no  indication  that  the  ratio  Xrms/f'b  is  ever  significantly  different  than  1 
for  this  shear  instability. 


1.  Introduction 

While  steaming  with  the  current  along  a  deep  tidal 
channel,  we  had  the  good  fortune  to  sample  a  kilometer- 
long  set  of  large  billows  (Plate  1).  As  the  billows  faded 
on  the  echo  sounder,  we  turned  the  ship  around  to  run 
upstream  back  through  them.  Because  of  the  increasing 
tided  currents  and  problems  retrieving  our  microstruc¬ 
ture  profiler,  we  fell  back  and  then  repositioned  several 
times.  But  the  current  swept  the  billows  downstream  by 
us,  allowing  us  to  collect  a  second  set  of  measurements. 

This  paper  is  a  detailed  presentation  of  our  fortuitous 
encounter.  Our  simultaneous  measurements  of  micro¬ 
structure,  acoustic  backscatter,  and  horizontal  currents 
provide  a  unique  realization  of  energetic,  naturally  oc¬ 
curring  billows.  The  turbulent  behavior  of  these  large 
instabilities  is  qualitatively  different  than  their  low  Rey¬ 
nolds  number  laboratory  counterparts;  specifically,  the 
billows  are  strongly  dissipative  from  very  early  in  their 
evolution  and  remain  actively  turbulent  for  up  to  6 
buoyancy  periods,  3  times  longer  than  found  in  labora¬ 
tory  studies  [Thorpe,  1973;  Koop  and  Browand,  1979]. 

Most  intriguing,  however,  is  the  change  in  properties 
along  the  kilometer-long  set  of  billows  that  must  result 
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from  changes  in  conditions  where  the  billows  first  form. 

The  mean  density  of  individual  billows  varies  within 
the  set;  we  propose  that  this  results  from  changes  in 
vertical  offset  of  the  velocity  and  density  interfaces.  The 
longitudinal  density  structure  we  observe  has  important 
implications  for  the  final  state  of  the  patch  ai  mixed 
fluid  produced  by  the  instability.  The  dissipation  rate 
within  the  billows  vuies  with  billow  density,  suggesting 
that  the  interface  offret  significantly  affects  the  intensity 
of  turbulence. 

Section  2  describes  our  instrumentation  and  then  the 
setting  and  spatial  organization  of  the  sampling.  We 
give  an  overview  of  the  data  in  section  3,  including  a 
simple  model  of  the  event  that  is  consistent  with  the 
observations.  Section  4  addresses  the  composition  of 
mixed  fluid  produced  when  the  shear  layer  becomes  tur¬ 
bulent  and  how  the  outer  scales  and  energy-containing 
scales  of  the  mixing  layer  vary.  We  summarize  our  re¬ 
sults  and  discuss  the  implications  in  section  5. 

2.  Background  cnALTTY  ■n'Tr^''”rrrHIT)  w 

2.1.  Instrumentation 

Profiles  of  velocity  and  acoustic  backscatter  were  col¬ 
lected  continuously  from  the  R/V  Miller,  a  15-m  utility 
boat.  A  150-kHz  RD  Instruments  acoustic  Doppler  cur¬ 
rent  profiler  (ADCP)  monitored  the  currents;  BioSonics 
120-kHz  and  200-kHz  narrow  beam  echo  sounders  inv- 
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Piute  1.  AioMsiic  l)a(  ks<  iitti  T  imaur  of  lli<-  Killows,  overlaid  with  profiler  1  rajeclories  (white 
lines),  will  ill  run  paralh'l  (o  ( he  sideloix-  r<  (urns  from  I  he  advanced  iiiicrosi  riirl  nre  proiiler  ( A  M  I’) 
(laint  red  lines)  I'lie  verlii  al  exa,i;^<'ral ion  is  .')  1 ;  pi.xel  si/e  is  0  KJ.')  in  in  the  vertical  and  every 
second  hori/ontally  ;\MI'  nnniherin^  f'ives  the  last  two  dii;its  of  drop  nninhers  referred  to  in 
the  ii  xt  Red  I  orrespoinls  to  stronger  hackscatter.  the  scale  being  volume  scattering  strength  in 
decihels  per  I  ;il’a  at  I  in  I'lie  (iistaiiee  scale  is  of  ship  motion  relative  to  I)  h  M)*a  (hi)  m)  The 
hillows  fade  with  time  as  they  oci  nr  in  decreasing  st  rat ilication 


aged  haclisi  alter  inleiisilv  The  free  f.il ling,  loo.sely 
tethered  advanced  ninrost  nn  I  lire  pro|il<  r  (AMP)  mea 
siiri'd  Vert  ical  jiri  dih  s  of  teiiiperat  nre,  conductivity,  ( i-n 
timi't<T  s  ale  \erlieal  slie.-ir  of  horizontal  ciirri’itls,  and 
cent  mnler  scale  vertical  temperat  nn-  gr.nlients  every 
■Jilt)  Itttl  III  along  the  ship  track 

Iliese  ihri-e  iiistriiiiieiit  systems  provide  distinctly 
ililferi-nt  vK-ws  of  the  shear  instahility  An  invalnahle 
nnaliased  two  dimensional  picture  of  the  water  (olimiii 
IS  olitailied  hy  the  echo  sounders,  without  which  wi- 
could  not  ideiilily  the  hillows  rile  ADCP  measures 
the  average  <  haii.gi-  in  the  velocity  held,  although  the 
averaging  tune  is  short  enough  that  the  large  coherent 
striiclnres  contrdnite  to  shorl  term  variat ions  in  the  vi'- 
locity  estimates  On  t  he  other  hand,  t  he  AMP  collects  a 
profile  of  |)oint  measnri'iiii'iits  (i c.,  a  small  volume  sam¬ 
ple  collected  in  a  fraction  of  a  .second)  relatively  cjtiickly. 
At  0  lif)  111  s  ',  the  proiiler  cro.s.ses  a  2()-m  billow  in 
about  dO  s,  a  lifti'eiitli  of  the  average  buoyancy  periotl 
across  the  first  large  billow  (A'  «  0.01-1).  AMP  profiles 
thus  strongly  ndlect  the  phase  of  tlie  coherent  struc¬ 
tures  they  sample,  we  examine  the  turbulence  within 
the  large  scale  structures  with  tiu'.se  data  ’Phis  variety 
of  perspei  lives  is  advantageous,  except  when  we  form 


variables  by  combining  fields  from  differi-nt  iiist rniiieiit 
systi'iiis,  owing  to  the  diffennce  in  averaging  limes, 
lh<-.s<-  variables  are  ipiile  “noi.sy,"  h-ading  to  a  problem 
akin  to  salinity  spiking  In  particular,  the  Hichardson 
number  is  formed  from  the  AMP  and  AIX’P  data  and 
appears  ipiile  noisy. 

I'ln-re  is  some  nncertaiiily  as  to  the  i-xact  ship  Ira 
ji'ctory  over  ground,  because  our  navigation  data  an 
poor,  'rriangnlatioii  from  radar  ranges  collected  ;is  th- 
.AMP  drops  Ix'gan  prfwides  navigational  lixi's  accurate 
to  about  50  III.  We  interpolate  between  known  |>osilions 
for  all  AMP  drops  lacking  fixes. 

rin-  ADCP  data  are  averaged  for  2  min.  producing 
ensembles  of  about  210  pings  I’he  averages  have  rms 
uncertainties  of  less  than  2  cm  s‘'  The  vertical  reso¬ 
lution  of  the  .\D(’I’  is  determined  by  the  1  in  bins  and 
8  m  pul.se  length,  which  approximate  a  Bartlett  filter 
having  a  d  dB  point  at  12  m  'Phe  vertical  shear  of 
the  horizontal  currents  is  formed  by  first-differencing 
the  -1  m  data;  =  (An/Ac)-  -p  (An/Ar)^,  where 
H  and  e  are  east  and  north  velocity  components,  re 
spectively,  and  Ac  is  -f  m  Phe  subscript  denotes  res¬ 
olution  of  scales  larger  than  12  m.  Unfortnnalely,  we 
are  unable  to  form  absolute  velocities  directly  from  the 
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ADCP  data  because  bottom-tracking  was  not  work¬ 
ing.  A  problem  with  the  compass  installation  on  the 
Aft//er  renders  the  measurements  unreliable  during  large 
changes  in  direction  and  noisy  at  other  times;  we  later 
note  when  this  is  a  problem.  We  use  a  one-dimensional 
linked-channel  tidal  model  to  predict  the  absolute  cur¬ 
rent  speed  [Lavelle  ei  ai,  1988].  Details  of  the  proce¬ 
dure  are  given  in  Appendix  A. 

The  AMP  was  dropped  from  the  stern  of  the  Miller 
while  under  way,  fell  at  0.6-0.7  m  s~^  to  within  roughly 
10  m  of  the  bottom,  and  was  retrieved  using  the  tether. 
Standard  processing  of  the  hydrographic  data  produces 
profiles  with  0.1-m  vertical  resolution.  Assuming  iso¬ 
tropy,  we  form  the  dissipation  rate  of  turbulent  kinetic 
energy  (e)  from  the  small-scale  shear,  and  the  dissipa¬ 
tion  rate  of  temperature  variance  (x)  from  the  vertical 
temperature  gradient  over  0.5-m  segments  as 

e  =  7.5i/  /  {2irk)^i^{k)dk  [W  kg-‘],  (la) 
Jko 

X  =  6kt  /  {2irk)^9T{k)dk  [K^  s-^j  (lb) 

Jka 

where  and  are  the  velocity  and  temperature 
spectra,  respectively;  v  is  the  kinematic  viscosity;  and 
ko  and  k^  are  the  lower  and  upper  bounds  of  the  integra¬ 
tion,  not  necessarily  the  same  for  the  two  integrals.  See 
Peters  et  ai  [1988]  and  Wesson  and  Gregg  [1994]  for 
details.  These  authors  estimate  that  absolute  errors  of 
ensemble  averages  are  roughly  a  factor  of  2.  Appendix 
B  is  a  brief  discussion  of  dissipation  rate  errors. 

Within  the  mixing  event,  much  of  the  water  column  is 
statically  unstable.  To  estimate  the  scale  of  overturns, 
we  reorder  the  density  profiles  to  be  monotonically  in¬ 
creasing,  a  procedure  first  proposed  by  Thorpe  [1977]. 
By  comparing  the  observed  and  reordered  profiles,  we 
form  profiles  of  displacement  (Lt),  that  is,  the  distance 
a  parcel  is  moved  within  a  profile  to  form  a  statically 
stable  water  column.  To  minimize  false  overturns  re¬ 
sulting  from  salinity  spiking,  we  sort  both  temperature 
emd  density  2uid  choose  the  sort  that  minimizes  Lt.  An 
overturn,  /lov,  is  the  distance  over  which  the  running 
sum  of  the  displacements  is  nonzero,  ^  Lt  #  0.  Many 
of  the  pareuneters  used  in  this  article  are  defined  as 


Lt 


Re  =  Auhfi/ 
Rcb  =  e/vN^ 
Rib  =  Apgh/p{AU)^ 
h  =  \AU/{du/dz) 

maxi 
Lt,  I 


^ov)  Xovi  Pm 


displacement  distance; 
overturn  height  (distance 
over  which  ELt  ^  0); 
layer  Reynolds  nund>er; 
buoyancy  Reynolds  number; 
bulk  Richudson  number; 
vorticity  thickness; 
thickness  of  elevated  e  layer, 
mean  e  over  L,; 
mean  overturn  properties; 

root-mean-square  displace¬ 
ment; 


Lb  =  buoyancy  (at  Ozmidov) 

scale. 

Profiles  of  JV*  are  formed  ftom  the  reordered  d^isity 
(p)  profiles,  ensuring  positive  values  of  N^.  We  regard 
these  estimates  to  be  most  representative  of  the  buoy¬ 
ancy  forces  against  which  the  turbulence  acts.  With¬ 
out  a  subscript,  is  formed  by  differencing  the  den¬ 
sity  over  1  m,  while  N^2  denotes  estimates  formed  after 
Bartlett  filtering  the  density  profiles  to  the  same  ver¬ 
tical  scales  as  the  ADCP  data  before  differencing  over 
4  m.  Estimates  of  the  gradient  Richardson  number, 
Rii2  =  N12/S12,  are  thus  formed  from  variables  with 
similar  bandwidths. 

Narrow  beam  echo  sounders  were  used,  having  a 
I2-m^  footprint  (circle  of  radius  2  m)  at  50-m  depth, 
much  smaller  than  the  50-  to  100-m  scale  of  the  billows. 
Profiles  were  collected  every  second  and  digitized  every 
1/8  m.  The  acoustic  system  is  calibrated,  allowing  us 
to  convert  images  to  volume  scattering  strength  (VSS) 
using  a  sonar  equation.  In  a  separate  study  we  com¬ 
pare  the  measured  VSS  with  estimates  of  backscatter 
from  the  small-scale  temperature  variance  measured  by 
AMP  (H.  £.  Seim  et  al.,  manuscript  submitted  to  the 
Journal  of  Atmospheric  and  Oceanic  Technology,  1993) 
(hereinafter  referred  to  as  Seim  et  al.,  1993).  The  es¬ 
timated  VSS  closely  matches  the  measured  VSS  at  the 
level  of  the  billows,  where  the  mean  temperature  gradi¬ 
ent  is  largest,  but  elsewhere  they  eue  uncorrelated.  We 
use  this  fading  to  justify  interpreting  the  upper  reflec¬ 
tive  horizon  at  this  depth  as  an  unaliased  image  of  the 
temperature  interface  (Plate  1). 

In  order  to  do  the  backscatter/microstructure  cor¬ 
relation  study,  it  was  necessary  to  estimate  the  posi¬ 
tion  of  AMP  on  the  echo  sounder  image.  Our  tech¬ 
nique  amounts  to  spatially  registering  the  AMP  profile 
on  the  backscatter  image  (Seim  et  al.,  1993),  as  shown 
in  Plate  1.  We  calibrate  the  model  with  the  appar¬ 
ent  echo  from  AMP  in  a  tremsducer  sidelobe  for  drops 
5855-5858,  which  appear  as  faint  red  lines  paralleling 
the  model  trajectories  in  Plate  1.  (Only  the  last  two 
digits  of  the  AMP  drop  numbers  are  given  on  the  plates 
emd  some  figures.)  We  estimate  that  calculated  trajec¬ 
tories  ue  within  ±15  s  of  the  true  trajectory  (see  Seim 
et  al.  (1993)  for  details). 

The  acoustic  images  are  collected  while  moving  and 
thus  represent  chamges  in  both  space  and  time;  they 
cannot  be  considered  a  “snapshot”  of  the  spatial  struc¬ 
ture  nor  a  time  series  of  the  evolution.  The  ship’s  mo¬ 
tion  distorts  the  horizontal  scales  by  Doppler  shifting; 
observed  features  of  length  L  have  a  true  wavelength  of 

A  =  L(1  -  U/Ump)  (2) 

where  the  flow  velocity  is  U-  Owing  to  shear  in  the  wa¬ 
ter  column,  the  Doppler  shift  varies  with  depth,  further 
distorting  the  shape  of  features. 
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2.2.  Setting 

The  data  were  collected  in  Admiralty  Inlet,  a 
30-lun-long  tidal  channel  con  :ng  Puget  Sound  and 
the  Strait  of  Juan  de  Fuca  (I  1).  The  latter  leads 
to  the  open  ocean.  Puget  d  receives  a  consid¬ 
erable  amount  of  runoff  which  drives  a  baroclinic  ex¬ 
change  flow  across  Admiralty  Inlet  and  makes  it  salt- 
stratified.  A  shear  layer  typically  exists  between  the 
inflowing  saline  bottom  layer  and  outflowing  surface 
layer,  providing  a  natural  laboratory  in  which  to  study 
shear  instability.  The  flow  is  complicated  by  the  2-  to 
4-m  mixed  (both  diurnal  and  semidiurnal)  tides  of  the 
northwest  coast  of  North  America,  which  in  combina¬ 
tion  with  the  exchange  flow  gives  rise  to  a  time-varying 
shear  flow. 

Our  samples  were  collected  off  Bush  Point,  a  lateral 
constriction  just  seaward  of  the  mouth  of  Hood  Canal. 


A  series  of  sections  was  collected  off  Bush  Point  dur¬ 
ing  and  after  an  extended  high  tide.  During  most  of 
these  sections,  wavelike  and  cat’s-eye-like  patterns  in 
the  water  column  appeared  on  the  echo  sounder.  We 
here  concentrate  on  the  sections  with  the  most  extensive 
microstructure  sampling  and  for  which  we  have  some 
navigation  data. 

The  underlying  cause  of  instabihty  at  this  location 
appears  to  be  linked  to  internal  hydraulics.  Bush  Point 
forms  a  pronounced  constriction  that  can  act  as  a  hy¬ 
draulic  control  [e.g.,  Armi,  1986].  The  influence  of  the 
topography  on  the  ebbing  flow  is  apparent  in  a  section- 
ally  averaged  tidal  transport  model  [Lavelle  et  ai,  1988] 
that  predicts  a  nearly  twofold  increase  in  current  speed 
as  the  flow  passes  by  Bush  Point.  Because  of  the  os¬ 
cillatory  tidal  flow,  however,  the  control  is  not  always 
active.  For  a  two-layered  flow,  criticality  is  determined 
by  the  composite  Froude  number. 


Longitude 

Figure  1.  Location  map  and  bathymetry  of  Admiralty  Inlet.  The  box  marks  the  study  area. 
Note  that  Bush  Point  ruid  Marrowstone  Island  constrict  northward  (ebbing)  flow,  making  this 
section  a  hydraulic  control  when  the  density  field  is  favorable  (see  text). 
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Figure  2.  Salinity  profiles  (a)  downstream  and  (b)  upstream  of  Bush  Point.  Numbering  shows 
the  order  in  which  the  profiles  were  collected.  The  persistent  interface  at  0.6  MPa  upstream  and 
large  changes  in  interface  height  downstream  are  consistent  with  critical  flow  conditions  at  Bush 
Point  during  the  sampling  period. 


g'hi  9'h2 


(3) 


with  Un  and  hn  the  speed  and  thickness  of  the  layers, 
respectively,  and  g'  =  g^p/p  the  reduced  gravity.  We 
observed  a  roughly  two-layered  flow,  the  surface  layer 
being  1  kg  m~^  lighter  and  flowing  0.5  m  s~^  faster 
than  the  lower  layer.  At  high  water,  when  there  is  no 
net  flow,  there  is  an  exchange  flow  (Ui  =  —U2  and  hi  = 
h2  =  50  m).  Then  =  0.25,  the  flow  is  subcritical, 
and  the  constriction  off  Bush  Point  does  not  act  as  a 
control.  As  the  flow  begins  to  ebb,  the  layer  velocities 
change  and  begins  to  increase.  We  find  that  the  flow 
will  be  critical  =  1)  at  Bush  Point  for  mean  flow 
speeds  {{Ui  +  U2)/2)  in  excess  of  0.43  m  s“'.  During 
our  observations,  the  flow  speed  increased  from  below 
this  value  to  well  above  it,  suggesting  we  observed  the 
tremsition  through  criticality  and  the  establishment  of 
hydraulic  control. 

Changes  in  height  of  the  density  interface  down¬ 
stream  of  Bush  Point,  but  not  upstream,  confirm  this 
interpretation.  In  a  set  of  profiles  collected  upstream  of 
Bush  Point,  before  and  after  the  billows  samples,  the 
interface  is  centered  around  0.6  MPa  (Figure  2b).  In 
profiles  collected  downstream  of  Bush  Point  the  inter¬ 
face  is  initially  at  0.6  MPa  (Figure  2a,  profile  3),  in¬ 
dicating  subcritical  conditions  throughout  the  region. 
The  interface  downstream  then  rises  from  0.6  MPa  to 
0.3  MPa  (profiles  3-7)  as  control  is  established.  The  bil¬ 


lows  were  observed  during  the  transition.  Samples  after 
the  billows  (profiles  8-11)  find  the  interface  still  raised 
downstream  of  Bush  Point  but  unchanged  upstream, 
confirming  hydraulic  control  of  the  internal  flow. 

During  the  transition  from  subcritical  to  supercritical 
flow,  the  velocity  and  density  interfaces  are  separated. 
The  density  interface  rises  first,  followed  by  the  velocity 
interface.  Because  shear  instabihty  occurs  at  the  level 
of  maximum  shear,  not  maximum  stratification,  the  dif¬ 
ference  in  vertical  location  of  the  two  interfaces  signifi¬ 
cantly  impacts  the  composition  of  the  billows  f<»med. 

Because  biUows  have  roughly  zero  phase  velocity  relar 
tive  to  the  mean  flow  [Ho  and  Hutm,  1984J,  an  instabil¬ 
ity  within  a  mean  flow  evolves  downstream  of  the  gen¬ 
eration  site  because  it  cannot  propagate  upstream.  If 
the  phase  speed  of  the  forcing  that  produces  Hi  <  0.25 
moves  downstream,  the  event  will  appear  as  a  transient 
at  a  fixed  spatial  coordinate.  In  the  x-z  plane,  the  bil¬ 
lows  will  all  be  at  a  similar  stage  of  development,  their 
growth  being  obvious  only  by  examining  changes  in  time 
(Figure  3a).  If,  however,  the  forcing  is  topogrs^hically 
fixed,  the  instability  grows  (mly  in  the  downstream  di¬ 
rection.  In  this  case,  x-z  plane  realizations  from  differ¬ 
ent  times  will  all  be  similar  (Figure  3b). 

Though  a  hydraulic  control  should  generate  a  spa¬ 
tially  growing  instability,  the  flow  acceleration  due  to 
the  tides  makes  the  flow  unsteady.  Having  never  re¬ 
sampled  the  Bush  Point  cross  secticm,  we  cannot  say 
whether  Hi  <  0.25  was  maintained  there  for  very  long. 


10,054 


SEIM  AND  GREGG;  OBSERVATIONS  OF  A  SHEAR  INSTABILITY 


z 


z 


Figure  3.  Sketch  of  two  types  of  spatially  growing 
instabilities.  The  plane  of  the  page  is  the  x-z  plane, 
with  time  increasing  into  the  page,  (a)  Transient  evo¬ 
lution;  billows  at  same  phase  of  development  at  all  x, 
evolving  only  with  time.  Dashed  lines  mark  position  of 
the  packet  as  it  is  advected  downstream,  (b)  Spatially 
growing;  billows  at  same  phase  of  development  for  all 
time,  evolving  only  with  x. 

We  have  insufficient  information  to  specify  the  form  of 
evolution  of  the  instability.  By  following  the  billows  at 
their  advection  speed,  though,  we  can  assess  their  evo¬ 
lution  in  time,  regardless  of  the  space- time  development 
of  the  shear  layer.  Not  knowing  if  the  instability  evolves 
spatially  or  temporally,  the  only  way  to  unambiguously 
measure  its  growth  or  decay  is  in  a  frame  of  reference 
moving  with  the  shear  layer.  We  adopt  this  coordinate 
system  when  assessing  evolution. 

Spatially  growing  instabilities  associated  with  con¬ 
strictions  were  found  by  Geyer  and  Smith  [1987]  in  the 
lower  reaches  of  the  Fraser  River  on  ebb  tide.  To¬ 
gether  these  studies  suggest  that  lateral  constrictions 
can  strongly  influence  vertical  mixing  in  channelized, 
stratified  flows. 

2.3.  Sampling  Pattern 

On  March  26,  1988,  soon  after  the  end  of  an  extended 
high  tide  (Figure  4a),  a  series  of  wavelike  features  was 
observed  at  middepth  on  the  echo  sounder  as  the  vessel 
steamed  northward,  with  the  ebbing  current,  past  Bush 
Point  (Figure  5).  The  ship  was  in  the  east  half  of  the 
channel  and  passed  over  a  small  rise  late  in  the  section, 
covering  3.4  km  over  ground  (Figure  4b).  We  designate 


the  set  of  AMP  profiles  collected  during  this  time  as 
section  A. 

As  the  acoustic  signature  of  the  features  faded,  the 
vessel  turned  southward  and  steamed  into  the  current. 
The  Miller  did  not  make  headway  and  was  carried 
north,  downstream,  while  steaming  into  the  current; 
the  four  profiles  collected  at  this  time  form  section  Al. 
After  repositioning  about  500  m  to  the  south,  three 
profiles  were  collected  while  slipping  northward,  paral¬ 
leling  section  Al;  these  we  designate  section  B.  Three 
more  profiles  were  subsequently  collected,  called  post-B 
(Figure  4b).  Table  1  summarizes  the  drop  numbers  and 
times  of  each  of  the  sections. 

The  ship  trajectory  relative  to  the  billows  is  found  by 
integrating  the  relative  velocity  between  the  ship  and 
the  water  at  0.5  MPa,  the  depth  at  which  the  billows 
are  centered  (Figure  4c).  We  estimate  the  ship  track 
during  the  few  minutes  between  drops  5858  and  5859 
when  the  Miller  reversed  direction  owing  to  a  problem 
with  the  compass  at  that  time  (marked  with  a  dashed 
line  on  Figure  4c).  Sections  Al  and  B  apparently  sam¬ 
ple  within  a  few  hundred  meters  of  the  fluid  in  which 
the  billows  were  found  in  section  A.  Section  Al  sam¬ 
ples  roughly  the  same  water  mass  at  the  north  end  ol 
section  A,  while  section  B  samples  near  the  first  bil¬ 
lows  seen  in  section  A.  The  later  drops  (post-B)  sample 
the  same  downstream  area  as  the  earlier  sections  (Fig¬ 
ure  4b)  but  well  after  the  billows  observed  in  section  A 
have  advected  downstream. 

In  summary,  the  Miller  steamed  3.4  km  downstreair 
relative  to  the  bottom  with  the  current  while  collecting 
section  A.  She  then  turned  into  the  ciurent  and  trier 
to  hold  position,  but  fell  off  while  profiling,  forcing  he; 
to  reposition.  Relative  to  the  flow  at  the  depth  of  th< 
billows,  the  Miller  sampled  a  1.8-km  segment  on  sec 
tiou  A,  then  turned  into  the  current,  such  that  the  fluir 
sampled  in  section  A  was  advected  past  the  ship.  I 
appears  that  sections  Al  and  B  sample  approximately 
the  same  fluid  as  in  section  A.  We  call  the  latter  refer 
ence  frame  billows-following  and  the  former  geographic 
They  are  surrogates  for  Lagrangian  and  Eulerian  refer 
ence  frames. 

3.  Observations 

This  overview  presents  the  observations  for  each  o 
the  sections  sequentially.  After  section  A  we  offer  ai 
interpretation  of  what  we  sampled. 

3.1.  Section  A 

3.1.1.  Before  billows.  As  section  A  begins,  th* 
water  column  is  continuously  stratified,  but  by  AMI 
5852,  two  moderately  stratified  layers  {N  »s  0.005  - 
0.015  rad  s"^)  are  separated  by  a  hi^-gradient  re 
gion  {N  «  0.03  rad  s“^)  (Figure  6).  The  mean  den 
sities  of  the  two  layers  differ  by  0.9  kg  m~^  with  abou 
a  0.5  kg  m“^  change  across  the  high-gradlcnt  region 
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Figure  4.  (a)  The  tides,  with  the  measurement  interval  shaded.  Sampling  begins  at  the  end 
of  an  extended  high  water  and  continues  through  maximum  ebb.  (b)  Enlarged  location  map, 
showing  the  eastern  half  of  the  channel  as  outlined  in  Figure  1.  Solid  symbols  mark  drops  whi^ 
passed  through  billows  or  large  overturns.  Triangles  mark  drops  from  section  A,  squares  mark 
section  Al,  circles  mark  section  B,  and  inverted  triangles  mark  post-B.  (c)  Progressive  vector 
displacements  of  ship  motion  relative  to  0.5  MPa,  wUch  presents  the  sampling  pattern  in  a 
billows-foUowing  reference  frame.  Dashed  line  marks  the  estimated  trajectory  discussed  in  the 
text.  Symbols  are  as  in  Figure  4b. 
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Figure  5.  Uncorrected  acoustic  image  of  sections  A-Al  showing  the  location  of  the  billows 
relative  to  bathymetry.  The  dark  region  at  the  base  of  the  image  is  the  bottom.  Note  cusping  of 
scattering  layers  above  the  billows  and  at  0.6  MPa  for  time  greater  than  1800  s.  We  interpret  it 
as  the  result  of  circulation  around  the  billows. 
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'Kible  1.  Section  Drop>s  and  Times  on  March  26,  1988 


Section 

AMP  Drops 

Times* 

A 

5849-5858 

2008-2105 

A1 

5859-5862 

2112-2130 

B 

5863-5865 

"141-2155 

post-B 

5866-5868 

2205-2232 

*Times  are  in  UTC. 

AMP,  advanced  microstucture  profiler. 


The  velocity  field  is  also  roughly  two-layered,  ejchibit- 
ing  weak  shear  in  the  two  layers  and  larger  shear  across 
the  high-density  gradient.  The  velocity  difference,  AU, 
between  the  layers  is  about  0.45  m  s“^.  Between  drops 
5849  and  5853,  the  gradient  region  between  the  two 
layers  thins  and  rises  from  below  0  6  MPa  to  above 
0.5  MPa,  during  which  time  doubles  to  8  x  10“^  s~^ 
and  Si2  triples  to  over  2  x  10~®  s"^.  The  echo  sounder 
reveals  a  nearly  horizontal  scattering  layer  at  the  depth 
of  the  high  region  which  also  thins  and  rises  (Plate  2). 

The  Reynolds  number  across  the  interface  is  very 
large,  Re  =  AUhft/  =  3.5  x  10®,  where  the  vorticity 
thickness,  h  =  \AU /{du I dz)max\,  is  11  m  in  AMP  5852 
and  1/  =  1.4  x  10“®  m^  s~^.  This  value  of  initial  Re  is 
comparable  to  the  Re  of  atmospheric  billows  observa¬ 
tions  and  is  100  times  larger  than  most  numerical  and 
laboratory  experimental  studies  of  stratified  shear  lay¬ 


ers  [Thorpe,  1987],  with  the  notable  exception  of  worl 
on  high-speed  shear  layers  [e.g.,  Brown  and  Roskko 
1974;  Dimotakis,  1989]. 

In  the  lower  portion  of  the  high  shear  layer,  Rin  < 
0.25  in  drops  5852-5853  (Figures  6  and  7),  owing  to  th< 
maximum  ATfj  layer  being  thinner  and  offset  upwarc 
from  the  layer  of  maximum  shear.  Setting  Ap  =  AU  = 
0.5,  the  bulk  Richaurdson  number  Rib  =  Apgh/p{AU)- 
rapidly  decreases  &om  0.38  AMP  5850  to  0.21  in  AMI 
5852  as  h  changes  from  20  m  to  11  m.  The  high-gradien 
region  therefore  satisfies  the  necessary  condition  Ri  < 
0.25  [Miles,  1961]  for  instability.  The  first  large-scalt 
billows  are  obvious  in  AMP  5854,  but  smaller  structure: 
are  visible  before  this,  beginning  at  AMP  5853. 

The  increase  in  Sjj  between  profiles  5850  and  585^ 
appears  to  “trigger”  the  billows.  The  10-m  rise  of  th« 
velocity  and  density  interfaces  while  passing  Bush  Poin' 
implicates  the  constriction  as  the  cause  of  the  shea 
enhancement.  A  thinning  interface  implies  dw/dz  <  0 
However,  in  the  contraiction  [y  downstream),  du/dx  i: 
also  negative,  requiring  a  large  downstream  acceleratioi 
{dv/dy  >  0)  to  satisfy  continuity.  Field  acceleration: 
may  provide  the  necessary  force. 

The  water  column  before  the  appearance  of  the  bil 
lows  is  far  from  quiescent.  At  all  depths  e  seldom  fall 
below  10“®  W  kg“^  and  typically  exceeds  10“®  W  kg“ 
at  several  depths.  The  buoyancy  Reynolds  number 
Reb  =  efvN^,  a  coiiunonly  used  measure  of  turbulen 
activity,  typically  exceeds  100  over  much  of  the  wate 


Figure  6.  AMP  5852,  collected  prior  to  the  appearance  of  the  large  billows:  (a)  density  (<r#), 
2-min  averaged  velocity  (u,v),  and  c;  (b)  (shaded)  and  5*2  (dashed);  (c)  Rii2,  with  the 
dashed  line  marking  0.25;  and  (d)  Reb  and  if  The  o-g  and  AT*  profiles  end  at  0.77  MPa  because 
of  a  clogged  conductivity  sensor  below  this  depth. 
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Plato  2.  Sunirnary  of  section  A.  All  panels  have  in  rominoii  the  bottom  axis  and  contours  of 
density.  Note  that  contours  from  AMI’  data  are  alia.sed  and  do  not  reflect  horizontal  feai tires 
<;I00  m  (a)  Acoustic  image  of  the  billows  overlaid  with  AMI’  trajectories  (b)  the  2  inin 
averaged  AD(’P  velocity  profiles.  Hlue  shading  shows  where  >  3  x  U)  *  s~",  and  yellow 
shading  shows  where  3  x  10“'^  <  >  10“'’  s“‘  (c)  Density  profiles  from  AMI’;  shading  is  for 

jVfT,  using  the  sanu'  levels  as  for  Sf,  (d)  Contours  of  log  \  (e)  I’rofiles  of  (f)  Contours  of 

log  (  In  Plate  2d  and  2f,  values  greater  than  -  7  are  shaded  with  a  contour  interval  of  1 
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Figure  7.  Profiles  (a)  as  features  appear  below  the  interface  (drop  5853)  and  (b)  through  first 
large  billow  (drop  5854).  The  various  panels  are  described  in  Figure  6,  except  vertical  velocity 
(w)  is  added  to  the  right-hand  panel  of  Figure  7b  (see  text). 

column,  well  above  the  minimum  value  of  30  needed  for  ei  al.,  1984].  At  the  level  where  is  largest,  iJet  <  10‘ 
the  turbulence  to  support  a  negative  buoyancy  flux  in  in  drops  5850-5852,  considerably  less  than  peak  values 
grid-generated  turbulence  [iZoAr  ei  al,  1988).  Where  greater  thM  10®  found  outside  the  high-gradient  re^oc 
iZcb  exceeds  200  the  turbulence  is  isotropic  at  dissipar  (Figure  6)  but  high  enough  to  support  active  mixing 
tion  scales,  and  where  Rcb  >  4  x  lO'*,  fully  developed  Displacement  profiles  indicate  the  high  gradient  regior 
isotropic  inertial  subranges  can  be  anticipated  [Gargett  remains  quite  stable  (Figure  6). 


SEIM  AND  GREGG:  OBSERVATIONS  OF  A  SHEAR  INSTABILITY 


10^ 


3.1.2.  Billows.  Begiimmg  at  AMP  5853,  the  mid¬ 
depth  scattering  layer  develops  wavelike  features  (Plate 
2)  which  quickly  grow,  appearing  to  overturn  as  AMP 
5854  is  collected.  We  interpret  the  features  as  billows 
because  of  their  location  at  a  sheared  density  interface. 
The  acoustic  signatures  of  the  billows  fade  as  the  ves¬ 
sel  moves  northward;  by  AMP  5858  they  are  no  longer 
visible.  As  discussed  below,  the  billows  occur  in  pro¬ 
gressively  less  stratified  fluid  to  the  north,  resulting  in 
rapidly  decreasing  values  of  x  (Plate  2),  consistent  with 
our  contention  that  the  backscatter  at  this  depth  re¬ 
sults  from  scattering  off  refractive  index  microstructure. 
Mimicking  where  is  largest,  by  AMP  5856  a  scat¬ 
tering  layer  develops  at  0.6  MPa  which  persists  through 
the  end  of  section  Al,  after  the  upper  scattering  layer 
fades  away. 

In  the  billows-following  frame  (Figure  4c)  the  aver¬ 
age  wavelength  over  the  first  four  bUlows  is  70  m.  If  the 
billow  advection  speed  equals  the  mean  shear  layer  ve¬ 
locity,  as  predicted  by  theory  and  observed  in  the  labo¬ 
ratory,  then  Doppler  shifts  should  be  removed  from  this 
estimate  of  the  wavelength.  The  peak-to-trough  ampli¬ 
tude  of  the  billow  sampled  by  AMP  5854  is  about  15  m 
on  the  echo  sounder  image,  increasing  to  a  maximum  of 
20  m  in  AMP  5855. 

Beginning  with  AMP  5854,  displacement  profiles  show 
a  large  overturning  region  (Plate  2c)  which  corresponds 
to  the  location  of  the  billows  on  the  echo  sounder. 
The  overturns  are  initially  below  the  density  interface 
(23.75  kg  m"®)  but  occur  above  it  by  the  end  of  section 
A,  as  demonstrated  by  overlaying  the  density  profiles  for 
AMP  5854-5858  (Figure  8).  Dissipation  profiles  show  e 
and  X  enhanced  over  the  same  5-  to  20-m-thick  sections 


Figure  8.  Overlaid  density  profiles  lor  AMP  5854- 
5858,  demonstrating  the  decrease  in  density  of  the  bil¬ 
lows  to  the  north  in  the  bUlows-following  frame. 


where  billows  appear  on  the  echo  sounder  and  the  di»- 
placemrat  profiles  reveal  large  overturns  (Plate  2).  The 
billows  are  also  apparent  in  profiles  of  x-  Even  more 
dramatically  than  e,  x  becomes  large  only  at  the  depth 
of  the  billows  (Plate  2d). 

The  velocity  interface  spreads  after  drop  5854,  caus- 
ing  Si2  to  decrease  (Plate  2b)  and  the  vorticity  thick¬ 
ness  to  increase  from  1 1  m  in  AMP  5852  to  22  m  in  AMP 
5858.  The  change  in  the  velocity  prcrfiles  is  consistent 
with  rapid  momentum  mixing  at  the  interface  beginning 
at  drop  5854.  It  would  be  surprising  ff  the  decrease  in 
shear  (thickening  of  the  shear  layer)  were  due  to  the 
increasing  channel  width.  Continuity  suggests  the  di¬ 
vergence  of  the  cross-stream  flow  would  be  balanced  by 
thinning  of  the  layer,  du/dz  =  —dv/dy  —  5w/dz  >  0 
or  dw/dz  <  0.  We  may  actually  underestimate  the 
vorticity  thickness  downstream  of  Bush  Point. 

The  acoustic  images  give  no  indication  of  vortex  pair¬ 
ing,  consistent  with  the  large  initial  Ri  as  found  by 
Koop  and  Browand  [1979],  although  Brown  and  Roshko 
[1974]  note  pedrings  are  difficult  to  identify  in  still  shots 
(they  use  cinefilms  to  identify  pairing  events).  The  data 
from  AMP  5854,  however,  strongly  suggest  the  profiler 
passed  through  two  billows  (Figure  7b).  The  vertical  ve¬ 
locity  within  the  billows,  estimated  from  perturbations 
in  AMP’s  fall  rate  (see  H.  E.  Seim  and  M.  C.  Gregg, 
manuscript  in  preparation,  1993)  (hereinafter  referred 
to  as  Seim  and  Gregg,  1993),  shows  two  maxima  and 
two  minima,  mimicking  the  <lispla.cement  profiles.  This 
is  inconsistent  with  a  single  vortex;  a  near-vertical  pro¬ 
file  through  a  spanwise  vortex  possesses  at  noost  one  sig¬ 
nificant  maximum  and  minimum.  The  sense  of  motion 
is  as  expected,  dense  fluid  rising  and  light  fluid  sinking. 
Interestingly,  simulations  by  Fritts  [1984]  of  pairing  at 
a  strong  density  interface  produce  a  density  field,  and 
billow  structure,  similar  to  that  measured  in  AMP  5854 
(Figure  9). 

To  summarize,  section  A  begins  in  a  roughly  two- 
layer  flow.  Passing  Bush  Point,  the  high-gradient  re¬ 
gions  of  both  the  velocity  and  density  fields  thin  and 
rise,  greatly  increasing  the  gradients  between  the  two 
layers  and  leading  to  Riu  <  0.25  in  the  lower  half  of 
the  shear  layer  in  AMP  5852  and  5853.  The  unaliased 
acoustic  image  demonstrates  the  interface  is  initially 
smooth  edong  the  flow.  Coherent  structures  appear  be¬ 
neath  the  interface,  where  Rii2  <  0.25  at  AMP  5853 
and  rapidly  grow  to  20-m  billows  by  AMP  5855.  AMP 
profiles  through  the  billows  reveal  large-scale  overturn¬ 
ing  of  the  density  field  within  the  billows  and  large  dis¬ 
sipation  rates.  As  we  move  north  through  the  set  of 
billows,  density  profiles  show  that  the  billows  occur  in 
progressively  less  dense  fluid.  This  lighter  fluid  is  also 
less  stratified.  Consequently,  turbulence  within  the  bil¬ 
lows  generates  less  x  to  the  north,  causing  the  acoustic 
signature  of  the  billows  to  fade  out.  Velocity  profiles 
indicate  that  shear  rapidly  decreases  and  the  vorticity 
thickness  increases  after  the  billows  form.  This  is  consis¬ 
tent  with  the  extraction  of  energy  from  the  mean  shear 
by  the  shear  instability  and  resulting  turbulence. 
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Figure  9.  Contours  of  density  Held  showing  a  pairing  event  at  a  strongly  stratified,  sheared 
interface  from  simulations  by  fhtts  [1984].  The  structure  is  similar  to  that  seen  m  the  acoustics, 
and  a  vertical  profile  through  the  paired  billows  in  Figure  9d  would  closely  resemble  AMP  5854. 


3.1.3.  Interpretation.  If  the  shear  layer  is  spa¬ 
tially  growing  and  statistically  steady  in  time  (in  a 
geographic  reference  frame),  the  density  of  the  billows 
would  rapidly  change  a®  they  move  downstream.  The 
billows  would  entrain  l£i>.ge  amounts  of  fluid  only  from 
the  upper  layer  and  would  grow  dramatically.  Because 
they  at  best  maintain  their  size  after  AMP  5855,  this 
model  of  the  evolution  is  not  physically  realistic.  Ap¬ 
parently,  conditions  are  cot  steady 

We  now  present  a  simple  cartoon  that  explains  the 
observed  longitudinal  density  structure.  The  billows 
form  as  critical  conditions  are  estabhshed  at  the  con¬ 
striction;  it  seems  reasonable  that  the  density  structure 
of  the  billows  is  related  to  this  time-dependent  process. 
Figure  10a  shows  the  two-layer,  inviscid  development  of 
controlled  flow,  based  on  a  model  of  a  similar  process 
over  a  sill  [Geyer,  1990].  Initially,  the  flow  is  subcritical 
throughout,  with  the  interface  at  height  zqI  the  inter¬ 
face  smoothly  shoals  at  the  constriction  (the  upper  layer 
is  assumed  to  be  active).  As  the  tidal  currents  acceler¬ 
ate  the  flow  to  critical  conditions  at  the  constriction,  a 
cusp  forms  that  subsequently  propagates  downstream, 
raising  the  interface  amd  establishing  a  growing  regime 
of  supercritical  flow  downstream  of  the  control. 

Now  consider  a  position  slightly  downstream  of  the 
constriction  (zq)  where  a  spatially  growing  shear  insta¬ 
bility  originates.  As  critical  conditions  at  the  constric¬ 
tion  are  established,  the  interface  at  xq  is  depressed 
relative  to  its  height  at  the  control,  causing  it  to  be 


centered  beneath  the  velocity  interface  (Figure  10b,  pro¬ 
file  1).  As  the  cusp  propagates  through  xq,  the  interface 
rises;  we  propose  that  this  occurs  more  rapidly  than  the 
velocity  field  can  respond,  so  the  density  interface  is 
temporarily  centered  above  the  velocity  interface  (Fig¬ 
ure  10b,  profile  3).  After  the  cusp  has  moved  far  down¬ 
stream  of  Xo,  the  two  interfaces  realign  at  a  higher  po¬ 
sition  in  the  water  column,  z\. 

Figure  10c  is  a  snapshot  of  the  density  field  down¬ 
stream  of  Xq  for  the  viscous  case  when  billows  form; 
At  time  1,  the  density  interface  is  beneath  the  velocity 
interface,  amd  billows  form  above  the  density  interface. 
At  time  2,  the  density  interface  has  risen  to  the  shear 
maximum,  and  the  billows  now  entrain  equad  amounts 
of  fluid  from  both  sides  of  the  density  interface.  By 
time  3,  the  density  interface  is  above  the  velocity  inter¬ 
face,  and  the  billows  form  beneath  it.  Fluid  from  time  1 
is  now  far  downstream  of  the  origin  (Figure  10c),  fluid 
from  successively  later  times  being  closer  to  the  ori¬ 
gin.  The  resulting  density  structure  is  in  good  agree¬ 
ment  with  our  observations  (Plate  2c)  and  resembles  a 
smoothed  version  of  Figure  10a  at  an  intermediate  stage 
of  development. 

Cusps  on  the  scattering  layer  near  the  end  of  section 
A  (Figure  5)  likely  result  from  billows  centered  above 
it  at  0.5  ^IPa.  This  is  consistent  with  our  mode!  that 
the  dens]  y  interface  was  centered  beneath  the  shear 
layer  when  these  billows  first  formed.  When  the  ve¬ 
locity  and  density  interfaces  eure  of  different  thicknesses 
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Figure  10.  (a)  Sketch  of  the  development  of  hydrauhc  control  at  a  constriction  as  the  tidal 
current  increases.  Initially,  the  flow  is  subcritical  (sohd  line),  the  interface  rising  only  at  the 
constriction.  As  critical  conditions  (G^  =  1)  occur  at  the  constriction,  a  cusp  forms,  which  sub¬ 
sequently  moves  downstreaun,  establishing  a  growing  region  of  supercritical  fiov  .  (b)  Conditions 
at  xo  where  the  shear  instability  originates.  The  density  interface  moves  from  below  the  velocity 
interface  as  the  cusp  first  forms  (1),  to  centered  at  the  same  location  (2),  to  slightly  above  (3). 
After  the  cusp  is  far  downstream,  both  interfaces  realign  at  a  new  position,  zj  (4).  (c)  The 
resulting  density  structure  downstream  of  xq.  Fluid  from  (1)  is  farthest  downstream. 

ind  oflset  a  fixed  amount,  recent  laboratory  studies  of 
bear  instabiUties  find  that  billows  form  at  the  level  of 
naximum  shear,  causing  the  density  interface  to  cusp 
or  Ri  «  0.25  [Lawrence  et  ai,  1991].  As  Ri  is  de- 
reased,  the  cusps  become  more  peaked  and  fluid  is  en- 
rained  across  the  interface  into  the  spanwise  vortices. 

Ve  speculate  that  changing  the  offset  between  the  in- 
erfaces  will  produce  similar  veuiations  in  the  amount 
>f  fluid  entrained  by  the  billows. 

i.2.  Sections  Al,  B,  and  Post-B 

The  longitudinal  density  structure  and  large-scale 
>verturning  of  the  billows  is  also  seen  in  the  later  pro- 
lies.  However,  because  ship  speed  varies  substantially, 


the  acoustic  imagery  is  difficult  to  interpret,  and  bil¬ 
lows  are  no  longer  obvious  on  the  echo  sounder.  We 
omit  these  images,  but  note  that  cusps  on  the  scattering 
layer  at  0.6  MPa  in  section  Al  (Figure  5)  and  the  per¬ 
sistence  of  structured  scattering  layers  at  0.4r-0.5  MPa 
in  section  B  suggest  the  continued  presence  of  coherent 
structures  at  iniddepth. 

In  section  Al,  large  overturns  occur  in  weakly  strat¬ 
ified  fluid  above  the  main  density  interface  at  23.75 
kg  m“®  (Figure  11).  Within  the  overturns,  c  is  compa¬ 
rable  to  values  at  the  end  of  section  A,  but  x  »  smaH, 
rising  above  10"^  K*  s“^  only  within  the  stratified  fluid 
beneath  the  overturns,  where  c  is  also  elevated. 

Section  B  samples  water  that  was  stratified  in  AMP 
5856  and  5855  of  section  A  (Figure  4c),  but  by  AMP 
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Figure  11.  AMP  5860  and  5862  from  section  Al.  Fields  in  each  pemel  are  as  in  Figure  4.  Large 
overturns  and  e  occur  in  weakly  stratified  fluid  above  the  main  density  interface;  high  f  persists 
through  the  interface  in  both  drops. 

5863  and  5864,  section  B  contains  an  obvious  20-tn-thick  Distinct  from  the  rest  of  the  profiles,  in  AMP  5865  the 
well-mixed  layer  whose  density  is  roughly  0.2  kg  m“®  density  field  is  stratified  over  much  of  the  region  of  ele- 
greater  than  in  section  Al.  The  large  overturn  is  cen-  vated  c,  possibly  the  result  of  collapse  of  the  biUows, 
tered  near  ae  =  23.75  kg  m~®  and  is  associated  with  although  c  still  exceeds  10~*  W  kg“^ 
elevated  c  (Figure  12a).  At  the  edges  of  the  well-mixed  Velocity  profiles  for  both  sections  continue  the  trend 
layer,  centered  at  0.5  MPa,  x  above  10~^  s"^.  of  weakening  shear  seen  in  section  A.  Both  5i2  and 
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Figure  12.  AMP  5863  and  5865  from  section  B.  Fields  in  each  panel  are  as  in  Figure  4.  Large 
overturns  are  centered  on  the  main  density  interface  in  AMP  5863.  In  AMP  5865  the  density 
field  is  more  stratified,  but  a  pronounced  e  layer  extends  over  28  m. 


Ni2  are  reduced  by  at  least  an  order  of  magnitude  at 
C.5  MPa  relative  to  conditions  before  the  billows  (Fig¬ 
ures  11  and  12).  The  Rit,  values  for  sections  A1  (0.4- 
0.6)  and  B  (0.2-0.5)  are  well  in  excess  of  0.33,  a  value 
that  Thorfe  [1973]  found  to  be  the  maximum  achieved 
in  his  laboratory  work  and  that  Gtytr  and  Smith  [1987] 
observed  to  be  the  preferred  value  during  ebb  in  the 


Fraser  River  estuary.  It  may  be  that  growth  cessation 
exhibits  some  Re  dependence,  and  shear  layers  with 
very  large  Re  persist  to  slightly  larger  Ri. 

The  three  additional  profiles  collected  after  section 
B  demonstrate  that  velocity  and  density  fields  have 
changed  considerably  (Figure  13).  Large  shear  is  now 
concentrated  above  0.4  MPa,  and  Rin  >  0.25  (although 
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Figure  13.  Profiles  5866-5868  collected  after  the  mixing  layer  is  no  longer  ^parent. 
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not  much)  in  this  depth  range  for  all  three  profiles. 
While  there  is  a  10-m  overturn  between  0.4  and  0.5  MPa 
in  AMP  5866,  the  other  two  profiles  do  not  display  any 
large  overturns  where  e  is  high  in  the  shear  layer.  Con¬ 
ditions  have  changed  enough  that  the  middepth  mixing 
layer  no  longer  exists. 

In  summary,  after  section  A  the  vessel  turned  into 
the  current,  roughly  staying  with  the  north  end  of  the 
set  of  billows  while  collecting  section  Al.  These  profiles 
reveal  large  overturns  centered  above  the  density  inter¬ 
face.  Repositioning  to  a  billows-following  position  close 
to  the  first  billows  in  section  A  finds  large  overturns 
centered  at  the  density  interface  in  section  B.  Both  sec¬ 
tions  are  consistent  with  the  density  structure  observed 
in  section  A.  Later  drops  at  the  same  geographic  posi¬ 
tion,  a  billows-following  position  upstream  of  the  billows 
in  section  A,  find  the  velocity  field  changed  and  no  large 
overturns,  indicating  the  mixing  layer  no  longer  exists. 

4.  Comparison  with  Laboratory  Studies 

The  initial  development  of  the  instability  is  similar 
to  that  found  in  the  laboratory.  With  the  minimum 
vorticity  thickness  of  11  m  observed  in  AMP  5852  be¬ 
fore  the  biUows  appear,  the  nondimensional  wavenum¬ 
ber  a  =  k{h/2)  =  (2ir/70)5.5  =  0.49.  This  is  some¬ 
what  higher  than  “typical”  atmospheric  observations  of 
a  rs  0.4  [Fritis  and  Rasiogi,  1985].  Our  value  of  a  is  also 
quite  close  to  the  most  unstable  mode  predicted  from 
linear  perturbation  theory  (0.42-0.49)  for  error  function 
and  hyperbolic  tangent  profiles  [Thorpe,  1971]. 

The  aspect  ratio  of  the  billows  is  0.14  using  a  bil¬ 
low  amplitude  of  10  m.  Thorpe  [1973]  found  the  aspect 
ratio  to  vary  with  initial  graidient  Richardson  number, 
from  0.6  for  Ri  =  0.05  to  0.2  for  Ri  =  0.15;  thus  our 
observations  are  consistent  with  an  initial  Ri  >  0.15. 


In  generatl,  then,  the  macroscopic  character  of  the  bil¬ 
lows  is  similar  to  laboratory  observations  and  previous 
atmospheric  observations. 

Turbulence  and  mixing  generated  by  the  billows  are 
less  likely  to  mimic  laboratory  flows.  Below,  we  first 
document  the  nature  of  mixing  by  the  instability  and 
confirm  the  similwty  to  high-speed  shear  flows  by  ex¬ 
amining  the  density  structure  within  the  shear  layer. 
We  then  show  that  the  billows  were  dissipative  as  soon 
as  they  appeared  on  the  echo  sounder,  but  that  shear 
spectra  at  dissipation  scales  take  some  time  to  develop 
a  —5/3  slope  in  the  inertial  subrange. 

We  also  see  if  the  event  is  well  represented  as  a  spa¬ 
tially  growing  instability.  To  facilitate  a  comparison,  we 
use  scales  common  to  laboratory  studies  of  shear  layers 
[Koop  and  Browand,  1979;  Dmiotakis,  1989].  The  com- 
puison  leads  us  to  speculate  about  the  importance  of 
an  offset  of  the  velocity  and  density  interfaces  in  deter¬ 
mining  the  magnitude  of  the  dissipation. 

4.1.  Density  Structure 

The  density  structure  of  the  shear  layer  provides  an 
important  clue  about  mixing  within  the  billows.  Heav¬ 
ily  averaged  profiles  smootUy  transition  between  the 
bounding  layers  [e.g..  Brown  and  Roahko,  1974],  im¬ 
plying  that  mixing  within  the  shear  layer  produces  a 
broad  range  of  mixing  products,  as  in  gradient  diflu- 
sion.  Instantaneous  data  collected  at  scales  ^proach- 
ing  the  Taylor  microscale  [e.g.,  Breidenihal,  1981]  are 
very  different,  finding  mured  fluid  of  only  one  composi¬ 
tion;  the  averaged  profile  shape  results  from  averaging 
in  the  presence  of  coherent  structures  which  occur  with 
decreasing  frequency  away  from  the  centerline.  This 
finding  prompted  a  phenomenological  model  <£  mixing 
in  high  Re  mixing  layers  [Broadwell  and  Breidenihal, 
1982]  which  differs  markedly  from  clasuc  gradimt  dif- 
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Figiire  14.  Probability  density  function  (pdf)  plots 
over  23.4  <  <  24.2  for  individual  drops,  positioned 

according  to  geographic  location.  Solid  squares  mark 
the  mean  density  of  overturns  with  hov  >  2  m  and  are 
seen  to  correspond  to  a  significant  mode  in  each  distri¬ 
bution. 

fusion  models.  We  next  test  which  model  our  data  best 
fit. 

The  probability  density  functions  (pdf)  of  the  fluid 
density  within  23.4  <  <t«  <  24.2  for  each  drop  are  shown 
in  Figure  14.  These  isopycnals  bound  the  fluid  layer  in 
which  the  shear  instabihty  occurs.  They  were  chosen 
because  they  are  local  minima  of  Lb  and  Lrms)  allowing 
us  to  define  a  fluid  layer  that  should  change  only  owing 
to  internal  processes.  AMP  5853,  collected  before  large 
billows  were  visible  on  the  echo  sounder,  is  a  bimodal 
distribution,  with  modes  at  23-54  and  24.15  kg  m”^ 
representing  the  bounding  layers.  In  AMP  5854,  a  third 
mode  develops  at  23.95  kg  m"^,  which  corresponds  to 
the  mean  density  of  the  overturn  through  which  AMP 
passes  (computed  separately),  marked  below  the  pdf 
with  a  solid  square  in  Figure  14.  For  all  subsequent 
drops  in  sections  A,  Al,  and  B,  the  mean  density  of  the 
overturns  (pov)  corresponds  to  a  significant  mode  in  the 
pdfs  which  lies  between  those  in  AMP  5853,  confirming 
the  similarity  to  high-iic  laboratory  flows. 

Unlike  the  model  of  Broadwell  and  Breidenthal  [1982], 
Pov  changes  moving  downstream.  The  overturn  mode 
occurs  at  progressively  lower  density  in  section  A,  is 
low  in  section  Al,  and  then  rebounds  in  section  B.  This 
is,  however,  consistent  with  our  model  of  the  density 
interface  rising  with  respect  to  the  velocity  interface 
over  time  at  the  site  where  the  flow  first  becomes  \m- 
stable.  This  suggests  that  an  offset  of  the  velocity  and 


density  interfaces  can  greatly  affect  the  crxnposition  of 
the  mixed  fluid  produced  in  a  shear  lay«.  The  chang¬ 
ing  initial  conditions  permanently  impact  the  develop¬ 
ment  of  the  shear  layer.  The  resulting  vertically  ho¬ 
mogeneous,  horizontally  varying  density  structure  has 
important  implications  for  the  final  state  of  the  shear 
layer.  We  pursue  this  in  the  discussion. 

4.2.  Turbulent  IVansition 

In  laboratory  flows,  the  distance  downstream  of  the 
splitter  plate  beyond  which  spatially  growing  shear  lay¬ 
ers  are  turbulent  is  called  the  small-scale  transition 
location  [Ho  and  Huerrt,  1984].  Beyond  this  point,  tur¬ 
bulence  is  fully  developed  and  molecular-scale  mixing 
proceeds  efficiently  [Dimotakis,  1989].  The  establish¬ 
ment  of  an  inertial  subrange  is  a  useful  criterion  for  as¬ 
sessing  if  transition  has  occurred.  In  shear  '  that 
are  initially  laminar,  the  slopes  of  shear  spe  mp- 

totically  approach  -1-1/3  from  large  negative  i,e.g., 
—2).  We  next  form  spectra  and  examine  their  shape  and 
“universality.”  Our  observations  are  some  of  the  first  to 
examine  transition  in  a  buoyancy-affected  shear  layer. 

Average  shear  spectra,  formed  over  10-  to  20-m  seg¬ 
ments  from  within  the  biUows,  approximate  the  em¬ 
pirical  Nasmyth  spectrum  [Oaiey,  1982]  of  the  aver¬ 
age  €  over  the  interval  and  the  theoretically  derivea 
spectrum  of  Panchev  and  Ktsich  [1969]  (Figure  15). 
(Although  the  oceanographic  community  has  adopted 
the  Nasmyth  spectrum  as  the  spectral  form  against 
which  to  compare  measured  spectra,  several  theoreti¬ 
cal  predictions  of  the  universal  spectral  form  exist.  Ap¬ 
pendix  C  presents  the  theoretical  spectrum,  a  compari¬ 
son  to  the  Nasmyth  spectrum,  and  that  of  Pao  [1965].) 
All  shear  spectra  from  AMP  5854  to  5865  are  broad- 
banded,  that  is,  possess  a  subrange  whose  slope  is  be¬ 
tween  0  and  -1-1/3  and  exhibit  viscous  rolloffs.  We 
cannot  resolve  the  low-wavenumber  portion  of  the  sub¬ 
range,  because  vehicle  motion  contaminates  the  shear 
measurements  at  scales  greater  than  2  m.  Wavenum¬ 
ber  spectra  normalized  by  the  Kolmogorov  velocity  and 
length  scales,  and  (c/»^)^^^,  are  remarkably 

similar. 

This  indicates  that  a  turbulent  subrange  exists  from 
the  time  we  first  observed  the  billows  on  the  echo 
sounder,  that  the  billows  axe  dissipative  from  early  in 
their  evolution,  and  that  they  are  not  destroyed  by  the 
turbulence.  This  bears  out  suspicions  voiced  by  both 
Thorpe  [1973]  and  Koop  and  Browcnd  [1979]  that  their 
results  on  high-iZi  flows  (greater  than  0.15)  may  not  ^ 
ply  to  energetic  oceanic  and  atmospheric  instabilities. 
They  found  the  shear  layer  at  large  Ri  to  be  a  short¬ 
lived,  inefficient  mixing  agent  which  quickly  collapsed 
under  gravity. 

In  our  case  the  instability  is  highly  dissipative  and 
long-lived.  Given  AMP  5855  and  5864  were  collected 
an  hour  apart  and  N  os  10~^  rad  s~^,  this  mixing  event 
generated  strong  turbulence  for  6  buoyancy  periods. 
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Figure  15.  Profiles  of  small-scale  shear,  temperature,  temperature  gradient,  and  shear  and 
dissipation  spectra  formed  over  the  depth  range  marked  on  the  profiles.  Dashed  lines  are  Panchev 
and  Kesich  spectra  for  comparison,  using  <  from  our  standard  processing  averaged  over  the  same 
depth  and  u  =  1.43  x  10~®  m^  s“^  Dotted  lines  are  the  Nasmyth  spectra,  (a)  Drop  5865; 
T  =  7.7  X  10-«  W  kg-‘,  iicb  =  1.9  x  10^  and  a  =  2.  (b)  Drop  5857:  ?  =  1.4  x  lO"®  W  kg-i, 
Rcb  =  9.2  X  10®,  and  a  =  1.6. 


This  is  3  times  longer  than  found  in  low-iZe  laboratory 
studies  [Thorpe,  1973;  Koop  and  Browand,  1979],  even 
for  much  lower  Ri. 

Two  departures  from  the  universal  spectrum  are  no¬ 
table.  AMP  5854  and  5855  have  spectral  slopes  of 
-1-1/6  over  1-10  cpm,  slightly  steeper  than  the  Nasmyth 
spectrum  (see  Figure  15a).  This  is  consistent  with  the 
asymptotic  approach  to  a  -1-1/3  subrange  observed  in 
laboratory  flows  [Ho  and  Huerre,  1984]  and  suggests 
the  turbulence  has  not  yet  etchieved  an  equilibrium 
state.  These  spectra  also  roll  off  at  lower  wavenum¬ 
bers  than  the  predicted  spectra;  for  the  Panchev  and 
Kesich  spectrum,  a  Kolmogorov  constant  of  2  rather 
than  the  typical  value  of  1.6  provides  a  better  fit  to 
the  measured  spectra.  The  other  oddity  is  less  eas¬ 
ily  understood.  Our  spectra  systematically  peak  at 
wavenumbers  half  the  peak  in  the  Nasmyth  spectrum; 
that  is,  if  the  Nasmyth  spectrum  is  given  by  ♦N(e.«'). 
our  data  are  better  fit  with  2$n(c/2,i^).  As  shown  in 
Appendix  C,  the  Panchev  ud  Kesich  spectrum  peaks  at 
a  lower  wavenumber  than  the  Nasmyth  spectrum.  We 
have  independently  confirmed  that  our  data  collection 
and  analysis  software  correctly  reproduces  an  electroni¬ 


cally  generated  white  noise  signal,  and  that  line  spectra 
are  properly  measured. 

In  AMP  5853  dissipation  rates  are  high  in  and  below 
the  density  step  at  0.45  MPa,  before  large  overturns 
form  (Figure  7a)  and  where  Riit  <  0.25.  Spectra  from 
this  drop  are  broadbanded,  but  the  rolloff  is  not  well  fit 
by  the  theoretical  spectra.  This  is  in  marked  contrast 
to  low-iJc  laboratory  experiments  in  which  turbulence 
was  observed  only  after  the  billows  attained  maximum 
amplitude  [Thorpe,  1973;  Koop  and  Browand,  1979]. 
However,  this  is  in  good  agreement  with  atmospheric 
observations,  in  particular  those  of  Browning  [1971].  He 
observed  17  “billow”  events  with  radar,  all  of  which 
were  preceded  by  a  turbulent  scattering  layer.  Indeed, 
radar  requires  the  presence  of  small-scale  turbulence  to 
image  the  luger-scale  instability  [Gage,  1989]. 

There  are  several  possible  reasons  for  the  discrepancy 
between  the  laboratory  observations  and  our  data.  In 
addition  to  the  luge  difference  in  /Ze,  an  effort  is  often 
made  in  laboratory  studies  to  keep  the  upstream  flow 
laminu,  whereas  in  the  tidal  channel,  the  flow  can  sel¬ 
dom  be  considered  laminu.  We  cannot  isolate  a  partic- 
ulu  aspect  of  the  flow  that  gives  rise  to  the  discrepancy. 
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but  note  that  tranaition  in  aome  geopbyaical  flowa  may 
be  qualitatively  different  from  tranaition  in  ideal  labo¬ 
ratory  aettinga. 

4.3.  Spatial  Evolution 

To  find  a  simple  description  of  the  evolution  of  the  in¬ 
stability,  we  compare  our  data  with  a  spatially  growing 
shear  layer.  Because  the  shear  layer  is  already  unstsr 
ble  by  AMP  5853,  we  assume  an  origin  upstream  of  it, 
at  AMP  5852.  We  normalize  the  downstream  distance 
and  layer  thickness  by  the  initial  vorticity  thickness, 
h  =  11  m.  To  quantify  the  dissipation  in  the  shear 
layer,  we  average  e  over  Lf,  roughly  the  distance  over 
which  f  >  10“^  W  kg”^  where  Sjj  is  large.  The  ele¬ 
vated  €  layer  is  consistently  bounded  by  several  meters 
of  fluid  in  which  e  is  well  below  the  threshold. 

The  growth  rate,  h/x,  is  roughly  constant  at  0(0.005) 
(dotted  line  in  Figure  16a),  suggesting  the  shear  layer  is 
slowly  but  s  teadily  growing.  This  growth  rate  is  a  tenth 
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Figure  16.  Growth  curves  assuming  a  spatially  grow¬ 
ing  instability.  The  origin  is  taken  to  be  AMP  5852, 
where  hi  fts  11  m.  (a)  The  h,  L,  (dashed  line),  and 
a  linear  fit  (dotted  Une).  The  trend  in  is  similar 
to  that  for  h  but  less  consistent,  as  expected  for  an 
unaveraged  measurement,  (b)  The  e  and  similarity  pre¬ 
diction  (dotted  line),  using  cq  =  8.8  x  10~^  W  kg~^  and 
7o  =  8-8  X  10“’^  W  kg“^  for  the  growth  rate  marked  in 
Figure  16a. 
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Figure  17.  Outer  scales  as  a  function  of  distance  north 
in  the  billows-following  frame  (cf.  Figure  4).  Different 
symbols  are  used  for  each  section:  crosses  for  section  A, 
pluses  for  section  Al,  and  open  circles  for  section  B. 

that  of  unstratified  flows  in  the  laboratory  [DtmotaJbs, 
1989]  but  similar  to  growth  curves  of  stratified  shear 
layers  for  Hi  >  0.1  [Tkorpt,  1973;  Koop  and  Bnwand, 
1979).  Changes  in  the  mean  velocity  field  are  consistent 
with  spatial  growth  at  fairly  large  initial  Ri. 

Changes  in  the  dissipation  rate  are  not  consistent 
with  a  spatially  growing  shear  layer.  In  a  self-similar 
shear  layer,  the  dissipation  rate  is  expected  to  be  in¬ 
versely  proportional  to  h,  such  that  edz  —  constant. 
The  doubling  of  h  between  AMP  5854  and  5858  should 
cause  7  to  decrease  by  a  factor  of  2.  After  peaking 
at  c  w  8.8  X  10“*  W  kg“*,  the  observed  downstream 
decrease  significantly  exceeds  that  expected  for  a  self- 
similar  flow,  shown  by  the  dotted  lines  in  Figure  16b. 
The  mean  dissipation  rate  near  x/hi  =  200  changes  dra¬ 
matically  between  sections  Al  and  B;  obviously  spatial 
growth  alone  cannot  e]q>lain  the  observations. 

As  already  mentioned,  the  shear  layer  was  not  steady, 
and  it  suggests  a  reference  frame  moving  with  the  bil¬ 
lows  would  be  more  appropriate  to  assess  the  evolution. 
Replotting  ?  versus  distance  in  the  billows-following 
frame  (Figure  4c)  finds  it  regularly  decreasing  north¬ 
ward  (Figure  17).  Changes  in  turbulence  levels  aloig 
section  A  appear  largely  related  to  the  changes  in  posi¬ 
tion  within  the  set  of  billows.  This  suggests  the  changes 
in  <  are  due  to  the  change  in  offset  of  the  velocity  and 
density  interfaces  along  the  set. 

In  the  billows-following  frame,  AMP  5863  resamples 
the  shear  layer  between  AMP  5855  and  5856.  For  these 
drops  the  change  in  ?  agrees  with  predictiou  for  a  self- 
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imilax  flow  (Figure  16b),  that  is,  both  points  lie  along 
the  predicted  curve.  That  changes  in  ?  of  fluid  tracked 
in  a  biUows-foUowing  frame  are  consistent  with  spatial 
growth  lends  support  to  our  suggestion  that  differences 
in  T  along  the  set  arise  from  differing  initial  conditions 
rather  than  varying  degrees  of  temporal  decay. 

To  crudely  estimate  the  impact  of  interface  offset  on 
the  dissipation  rate,  we  assume  all  billows  have  the  same 
rate  of  downstream  decrease  in  e  (parallel  dotted  Unes  in 
Figure  I6b).  If  the  origin  is  the  same  for  aU  the  billows, 
the  initial  value  of  7  must  have  varied  by  slightly  more 
than  a  factor  of  10.  All  other  initial  conditions  being  the 
same,  which  of  course  is  a  big  assumption,  the  interface 
offset  may  lead  to  1  order  of  magnitude  variation  in  the 
maximum  value  of  7. 

5.  Overturns  Within  the  Shear  Layer 

We  briefly  characterize  turbulence  within  the  shear 
layer  using  scales  and  parameters  common  to  oceano¬ 
graphic  studies.  We  follow  Wijesekera  ti  ai  [1993]  and 
H.  Peters  et  al.  (The  mechanisms  of  equatorial  turbu¬ 
lent  mixing,  2,  Detail  and  scaling  of  mixing  events,  sub¬ 
mitted  to  the  Journal  of  Geophysical  Research,  1993) 
by  conditionally  sampling  only  overturns.  Within  the 
overturns,  values  aie  averaged  before  forming  the  var¬ 
ious  scales,  the  idea  being  that  turbulent  length  scales 
are  properly  defined  only  within  the  unstable  portions 
of  the  water  column. 

Every  AMP  drop  collected  after  the  billows  form 
shows  large  overturns,  hoy  >  5  m.  Of  our  12  profiles, 
none  sample  a  stable,  high-gradient  region.  Either  the 
profiler  never  passed  through  braids,  or  the  braids  in 
this  realization  are  also  sites  of  overturning.  Given  the 
large  initial  Re,  secondary  Kelvin-Helmholtz  instabili¬ 
ties  are  anticipated  [Corcos  and  Sherman,  1976],  as  are 
streamwise  vortices  [Breidenthal,  1981].  Evidence  for 
secondary  K-H  instability  can  be  seen  in  Plate  1,  where 
a  small  billow  is  imbedded  in  the  scattering  layer  just 
upstream  of  the  AMP  5854  trajectory. 

We  form  average  properties  for  all  overturns  greater 
than  2  m;  smaller  overturns  contain  too  few  data  points 
to  form  stable  estimates.  Considering  the  largest  over¬ 


turns  in  each  profile,  only  <ov  »nd  Noy  display  mono¬ 
tonic  changes  in  the  billows-following  reference  frame. 
Both  are  smallest  to  the  north  (section  Al);  the  trend 
in  (ov  mimics  that  of  7  from  the  previous  section.  Low 
N^„  in  section  Al  probably  reflects  that  most  of  the 
overturns  occur  in  weakly  stratified  fluid. 

Table  2  presents  the  mean  and  range  of  properties  of 
2dl  the  overturns  in  the  shear  layer.  For  all  the  observed 
overturns,  Re^  >  2000,  easily  satisfying  the  criterion  of 
Gargett  ei  al.  [1984]  for  isotropy  at  dissipation  scales. 
The  large  variations  in  poy  reflect  chuges  in  the  vertical 
offset  of  the  density  and  velocity  interfaces  over  the  set 
of  billows.  With  the  exception  of  the  largest  overturn 
in  drop  5861,  f  >  10“^  W  kg“*.  The  smaller  overturns 
are  quite  energetic  and  for  AMP  5858  and  5861  are 
associated  with  higher  Cov  than  the  larger  overturns. 
The  mean  overturn  size  within  the  shear  layer  is  7.1  m 
with  Lrm*  =  2.58  and  Lb  =  2.71.  The  ratio  Lrm»/^  = 
0.95  is  close  to  1,  as  has  been  observed  in  almost  all 
previous  oceanic  thermocline  studies  [e.g.,  Dillon,  1982; 
Peters  et  al.,  1988]. 

The  equivalence  of  Lb  and  Lnn>  is  often  taken  to  indi¬ 
cate  that  stratification  has  limited  the  maximum  scale 
of  overturning  [e  g.,  Rohr  ei  al.,  1988],  a  state  that  Ivey 
and  Imberger  [1991]  find  to  be  maximally  efficient  at 
mixing.  Using  the  largest  overturns  in  each  profile, 
we  find  no  regular  trend  in  Lb/Lrm*  %  1  with  time 
or  in  the  billows-following  frame,  although  Lmt  > 
in  section  A,  while  Lb  >  Lmt  in  section  B.  A  plot 
of  Lrm*  versus  Lb  reveals  wide  scatter  (Figure  18), 
0.5  ^  Lb/Lfnit  ^  4.3.  For  hoy  ^  5  m.  Lb  tends  to 
exceed  Lnn»;  for  hoy  >  5  m,  the  reverse  is  true. 

The  data  allow  estimation  of  the  mixing  efficiency, 
7  =  Xpe/cov .  where  Xpe  is  the  dissipation  rate  of  poten¬ 
tial  energy.  Because  Xpe  ceuinot  as  yet  be  measured  in 
a  salt-stratified  fluid,  we  estimate  it  from  Xi  assuming 
that  the  dissipation  rates  of  heat  and  salt  are  in  pro¬ 
portion  to  their  respective  contributions  to  the  meeui 
density  gradient.  Following  Gregg  [1987],  we  use 

Xpe  =  [W^]*(l  +  l//??)X  [Wkg-1]  (4) 

where  a  is  the  thermal  expansion  coefficient,  R^  = 
a{ffr/dz)/0(dS/dz)  ta  -0.05  is  the  stability  ratio,  T 


Table  2.  Average  Properties  of  the  25  Overturns  with  ^ov  ^  2  in 


Variable 

Mean 

Range 

m 

7.1 

2-19 

Pe,,  kg 

23.7 

23.5-24.13 

'^rniJr  in 

2.58 

0.43-7.45 

Lb,  nn 

2.71 

0.54-8.7 

Cov.  W  kg-‘ 

1.83  X  10-* 

5.9  X  10-*-8.4  X  10~* 

N\  s-» 

6.7  X  10~® 

4.8  X  10“*-3.1  X  10~* 

Rth,  negative 

3.0  X  10* 

2.2  X  10®-1.5  X  10® 

X„,  K*  8“' 

1.4  X  10“^ 

2.3  X  10“*®-1.2  X  10"* 

Xp.,  w  kg-‘ 

9.8  X  10"^ 

3.7  X  10"*-6.2  X  10~* 

7,  negative 

0.58 

0.17-1.3 
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Figure  18.  Scatter  plot  of  inni  versus  Lb  for  overturns 
greater  than  2  m  within  the  mixing  layer.  Overturns 
greater  than  5  m  are  marked  by  crosses,  those  with 
/lov  <  5  m  are  marked  by  o’s.  A  least  squares  fit  yields 
Li  =  0.86Lrmt>  vrith  a  correlation  coefficient  of  0.47,  or 
Kendall’s  r  of  0.49. 

is  temperature,  5  is  salinity,  and  0  is  the  haline  con¬ 
traction  coefficient.  For  the  largest  overturns,  0.17  < 
7  <  1.3,  with  an  average  of  0.58,  which  is  larger  than 
previous  estimates  of  7  based  on  simultaneous  €  and 
X  measurements  (0.24  [Oakey,  1982];  0.12  [Peters  and 
Gregg,  1988];  and  0.12  <  7  <  0.48  [Movm  ei  al,  1989]). 
Averaged  over  sections,  7  is  greatest  in  section  A  (0.72), 
least  in  A1  (0.35),  and  intermediate  in  B  (0.65).  This 
suggests  that  offset  of  the  interfaces  also  affects  7,  but 
uncertainties  in  Xpe  larger  than  the  differences. 

6.  Summary  and  Discussion 

6.1.  Summary 

We  made  two  passes  through  a  kilometer-long  set  of 
20-m-tall  billows  within  an  unstable  shear  layer  over 
a  time  period  of  1  hour.  Continuous  velocity  profiles 
indicate  a  rapid  increase  in  shear  passing  by  Bush  P^t. 
Initial  conditions  of  the  mixing  event  are  (1)  AU/2U  » 
0.5,  Aplflo  »  5  X  10~‘*  across  the  high-gradient  region 
in  a  continuously  salt-stratified  (large  Schmidt  number) 
flow.  The  initial  vorticity  thickness  is  A,-  »  11  m.  (2) 
Ri  <  0.25  just  beneath  the  large  density  gradient  before 
the  echo  sounder  reveals  billows.  (3)  »  3  x  10*. 

The  wavenumber  and  aspect  ratio  of  the  bUlows  which 
form  are  consistent  with  the  initial  observed  h  and  Ri. 
Within  the  mixing  layer,  several  more  points  can  be 
made. 

1.  Instantaneous  density  profiles  through  the  bil¬ 
lows  reveal  large-scale  overturning,  and  we  find  these 


overturns,  in  all  cases,  to  be  associated  wiUt  elevated 
dissipation  rates.  In  the  billowa-foUowing  liame,  the 
overturns  occur  in  progreanvely  less  dense  fluid  to  the 
north.  This  implies  changes  in  the  offMt  of  the  velocity 
and  density  interfaces  over  time  where  the  billows  first 
form. 

2.  The  instantaneous  distribution  of  density  is  con¬ 
sistent  with  mixing  models  of  turbulent  shear  layers  at 
large  Re  [Eroadwell  and  Breidenthal,  1982].  The  persis¬ 
tence  of  overturn  modes  throughout  the  mixmg  event, 
long  after  the  initial  formation  of  the  billow  and  de¬ 
spite  an  obvious  longitudinal  structure  in  the  density 
field,  supports  the  notion  that  billows  permanently  re¬ 
tain  characteristics  inqiarted  at  their  formation,  even  in 
a  continuously  stratified  fluid. 

3.  The  transition  to  turbulence,  its  persistence,  and 
the  mixing  efficiency  of  billows  in  geophysical  flows  can 
be  qualitatively  different  than  in  the  laboratory.  The 
billows  are  strongly  turbulent  over  the  entire  sampling 
period,  even  in  their  early  development.  Turbulence 
within  the  shear  layer  is  long-lived  in  compavison  to 
low-i2e  laboratory  experiments  which  find  th*  shear 
layer  to  be  relaminarized  within  two  buoyancy  periods 
[Thorpe,  1973;  Koop,  1976].  Our  observations  indicate 
that  strong  dissipation  persists  for  more  than  six  buoy¬ 
ancy  periods. 

4.  Although  h  increases  regularly  over  sections 
A-B,  like  a  spatially  growing  stratified  shear  layer,  I 
decreases  more  rapidly  than  1/h  in  section  A  and  in 
section  A1  is  much  less  than  expected.  Values  of?  from 
section  B,  which  resample  part  of  section  A,  roughly  fol¬ 
low  values  predicted  for  a  self-similar  shear  layer.  We 
conclude  that  spatial  growth  adequately  represents  the 
outer  scale  variations  but  that  changes  in  interface  off¬ 
set  strongly  affect  the  level  of  dissipation. 

5.  The  energy  scales  Lb  and  Lnm  display  no  sig¬ 
nificant  trends  in  geographic  coordinates  or  a  billows- 
following  reference  frame.  Their  ratio  remains  within 
a  factor  of  2  of  unity  for  the  largest  overturns  and  av¬ 
erages  0.95.  Turbulent  activity  of  overturns  remains 
strong  over  the  entire  sampling  period,  Ret  >  2000, 
and  large  billows  may  be  quite  efficient  mixers,  7  »  0.6. 

6.2.  Discussion 

Previous  studies  of  buoyancy-affected  shear  layers 
[Thorpe,  1973;  Koop  and  Browand,  1979]  have  empha¬ 
sised  the  inqiortance  of  Ri  in  determining  the  character 
of  a  shear  instability.  This  study  indicates  the  impor¬ 
tance  of  Re  and  interface  offset.  As  noted  by  Koop  and 
Browand,  the  parameter  space  to  be  explored  is  indeed 
large. 

Observations  by  Browning  [1971]  suggest  interface 
of&et  is  a  conunon  feature.  In  most  of  the  17  atmo¬ 
spheric  billow  events  he  observed,  the  interfaces  were 
oflset;  eight  billosrs  occurred  "on  the  edge  of  a  stability 
maximum”  and  four  occurred  in  uniform  stratification, 
although  all  were  observed  at  a  level  of  m*Timiiin  shear. 
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Together  with  our  data,  this  implies  that  displaced  ve¬ 
locity  and  density  interfaces  at  kilometer  scales  (the  size 
of  the  mixing  patch)  are  the  rule  rather  than  the  excep¬ 
tion  in  geophysical  fluids. 

The  long-lived  nature  of  active  turbulence  in  this 
study  implies  a  fundamental  difference  between  this 
natural  shear  layer  and  those  at  low  Re  in  the  labora¬ 
tory.  Energetic,  naturally  occurring  instabilities,  which 
transition  to  turbulent  motion  soon  after  their  forma¬ 
tion,  will  necessarily  be  more  effective  at  producing 
mixed  fluid,  especially  if  they  are  longer  lived  than  their 
small  counterparts.  The  mixing  efficiency  of  big  billows 
may  be  substantially  greater  than  that  of  small  billows. 
There  is  mounting  evidence  that  mixing  rates  are  very 
small  in  the  bulk  of  the  ocean  thermocline.  Our  study 
implies  that  the  importance  of  sites  of  vigorous  mix¬ 
ing  may  be  underestimated  if  only  their  scale  is  con¬ 
sidered;  they  may  also  be  more  efficient  at  producing 
mixed  fluid. 

The  first-order  variability  in  dissipation  rate  of  the 
shear  layer  is  captured  in  the  billows-following  refer¬ 
ence  frame.  However,  the  change  in  turbulence  levels 
between  sections  A  and  AMP  5864  and  5865  is  signifi¬ 
cant  and  suggests  the  turbulence  is  beginning  to  decay. 
If  this  is  the  case,  then  some  of  the  trend  in  e  with  layer 
density  may  be  due  to  different  sampling  times.  In  the 
billows-following  frame,  we  first  sampled  the  southern 
portion,  when  these  billows  were  first  forming.  We  sam¬ 
pled  the  northern  part  well  after  they  formed.  If  the 
billows  did  not  behave  in  a  self-similar  fashion,  we  have 
misinterpreted  the  data. 

We  did  not  track  the  mixing  event  long  enough  to 
observe  the  final  state  of  the  mean  fields.  In  this  study, 
all  but  the  last  profile  through  the  shear  layer  reveal 
a  large  mixed  layer  bounded  by  high-gradient  regions, 
similar  to  the  model  of  Woods  and  Wiley  [1972].  How¬ 
ever,  two  points  cloud  this  finding.  First,  in  AMP  5865 
the  density  field  appears  to  be  relaxing  to  a  single,  thick 
high-gradient  layer.  It  is  possible  that  up  to  this  point, 
coherent  structures  dominate  the  density  structure,  but 
that  after  they  are  dissipated,  unmixed  fluid  invades 
the  biUow  structures  and  the  density  field  relaxes  to  a 
smoothly  varying,  stable  state.  Second,  even  if  a  mixed 
layer  persists  after  the  billows  dissipate,  the  longitudi¬ 
nal  density  structure  in  the  shear  layer  will  generate  a 
circulation  that  will  tend  to  restratify  the  shear  layer. 
A  vertical,  homogeneous,  horizontally  stratified  layer 
will  generate  a  lock-exchange  type  of  flow,  tending  to 
produce  a  smoothly  varying  vertical  density  gradient. 
Further  observations  are  required  to  resolve  this  issue. 

If  we  estimate  e  as  q^/l,  taking  q  »  0.1A[/  and  /  = 
20  m,  and  use  the  mean  stratification  =  10“^  s”*, 
we  find  Rcb  =  4.5  x  10'*,  close  to  the  mean  value  of 
the  largest  overturns.  Apparently  simple  scaling  argu¬ 
ments  provide  good  order-of-magnitude  estimates  of  the 
turbulent  activity  of  the  mixing  event  regardless  of  the 
interface  offset. 


Unanswered  is  a  quantification  of  the  import  wee  of 
turbulence  in  an  energy  balance  of  the  mixing  event.  In 
particular,  we  can  assess  the  importance  of  dissipation 
in  changing  the  shear  and  how  much  of  the  change  in 
potential  energy  is  due  to  the  buoyancy  flux.  The  en¬ 
ergy  fluxes  into  and  out  of  the  shear  layer  provide  clues 
about  the  generation  mechanism  of  the  instability  and 
a  measure  of  energy  loss  due  to  internal  wave  radiation. 
These  topics  are  pursued  by  Seim  and  Gregg  (1993) 

Appendix  A:  Barotropic  Velocities 

Owing  to  poor  navigation  data,  we  are  unable  to  form 
absolute  velocities  with  the  ADCP  data.  We  estimate 
the  flow  speed  during  the  time  we  were  sampling  with 
a  linked  tidal  transport  model  of  Puget  Sound  devel¬ 
oped  by  Lavelle  et  al.  [1988].  Using  finite  difference 
techniques,  they  solve  a  set  of  linearized  momentum 
and  continuity  equations  in  one  dimension,  constrained 
to  be  the  best  fit  to  observed  water  level  data.  The 
model  predicts  the  amplitude  and  phase  of  the  princi¬ 
pal  tidal  components  at  over  50  calibration  points  to 
within  ±2  cm  and  ±5‘.  The  model  also  predicts  the 
average  velocity  at  589  cross  sections  throughout  the 
Sound.  While  this  should  be  an  equally  accurate  esti¬ 
mate  of  the  sectionally  averaged  velocity,  the  few  sec¬ 
tions  that  have  been  well  instrumented  display  large 
vertical  auid  horizontal  gradients  in  velocity  at  tidal  fre¬ 
quencies.  The  model  velocities  are  therefore  of  limited 
accuracy  at  a  given  lateral  position.  Interestingly,  Bush 
Point  is  one  of  the  few  sections  that  has  been  instru¬ 
mented.  Three  moorings,  each  with  multiple  current 
meters,  were  mauntained  across  Bush  Point  for  several 
months  in  1977.  In  cotidal  charts  of  the  velocity,  the 
shear  at  Bush  Point  is  principally  in  the  vertical  [Lavelle 
et  al,  1988].  The  secondary  lateral  structure  displays  a 
velocity  core  centered  in  the  east  half  of  the  channel. 

We  estimate  the  current  speed  at  our  AMP  drop  lo¬ 
cations  by  first  predicting  the  tidal  current  at  all  model 
points  between  the  sections  that  bound  the  area  we 
sampled  (Figure  Al).  We  then  interpolate  between  the 
curves  using  our  time  and  position  information.  Only  in 
section  A  is  the  interpolated  curve  much  different  than 
the  individual  curves.  For  this  section  the  flow  rapidly 
accelerates  as  it  moves  through  the  lateral  constriction. 
Using  data  collected  on  anchor  just  north  of  Bush  Point 
in  May  1988  with  the  same  instrument  system,  we  find 
the  model  underpredicts  the  speed  by  20%. 

Appendix  B;  Dissipation  Rate  Errors 

Assessing  the  uncertainty  associated  with  dissipation 
rate  estimates  is  not  straightforward.  Because  of  the 
highly  skewed  nature  of  the  distribution  of  e,  standard 
statistical  procedures  that  assume  normal  distributions 
cannot  be  applied.  Briefly,  we  consider  two  forms  of 
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Figure  Al.  Predicted  tidal  currents  at  the  appropriate  model  sections,  labeled  a-d,  moving 
from  south  to  north  through  Bush  Point.  The  solid  line  is  the  interpolated  current  along  the  ship 
track,  with  AMP  drop  numbers  marked. 


error  in  our  estimates:  a  bias  in  absolute  estimates  due 
to  calibration  and  processing  errors  and  a  random  er¬ 
ror  due  to  natural  variability  of  the  dissipation  which 
controls  the  precision  of  our  estimates.  We  are  less 
concerned  with  the  former  in  this  study,  because  the 
two  shear  probes  in  AMP  were  not  changed  during  the 
course  of  the  experiment  and  thus  a  calibration  error 
affects  all  data  equally.  As  stated  in  the  text,  we  be¬ 
lieve  we  are  within  a  factor  of  2  of  the  ^K;tu^d  dissipation 
rate. 

Of  more  concern  is  to  realistically  bound  estimates 
developed  for  each  AMP  profile  to  determine  whether 
the  trends  we  observe  are  statistically  significant.  This 
may  be  approached  in  a  number  of  ways.  F>om  theory, 
we  anticipate  that  e  is  lognormally  distributed  in  sta¬ 
tistically  homogeneous,  isotropic  turbulence  [e.g.,  Gut- 
vich  and  Yaglom,  1967].  Confidence  limits  for  lognor¬ 
mally  distributed  data  are  a  function  of  the  number 
of  samples  and  the  standard  deviation  of  the  natural 
logarithm  of  the  e  estimates,  o'Inc  [Land,  1975].  Using 
the  0.5-m  estimates  from  within  the  elevated  e  region 
of  the  shear  layer  (£.,),  we  find  o'Inc  <  1-5.  Over  20  m 
we  obtain  eighty  0.5-m  estimates  by  overlapping  data 
windows  50%,  which  suggests  the  confidence  limits  ue 
somewhat  less  than  a  factor  of  2.  The  validity  of  this  ap¬ 
proach  is  questionable,  since  many  of  the  distributions 
depart  from  lognormality  and  the  individual  estimates 
are  certainly  not  independent  (i.e.,  uncorrelated).  In  a 
practical  sense,  the  e  estimates  are  based  on  summing 
a  large  number  of  spectral  estimates,  each  of  which  is 
a  sum  of  squares,  and  hence  are  Xn  distributed  with  n 
degrees  of  freedom.  If  we  are  grossly  conservative  and 
assume  each  0.5-m  estimate  has  2  degrees  of  freedom,  a 


20-m  average  has  80  degrees  of  freedom,  or  confidence 
limits  of  about  20%.  The  precision  of  our  estimates 
likely  lies  somewhere  between  these  sets  of  limits. 


Appendix  C:  Shear  Spectra  at 
Dissipation  Scales 


Several  theoretical  predictions  of  the  universal  spec¬ 
tral  form  exist.  Pao  [1965]  presents  a  simple  model  of 
the  spectrum,  based  on  the  assun^tion  that  the  trans¬ 
fer  of  energy  au:ro6s  a  given  wavenumber  k  in  the  iner¬ 
tial  subrange  is  equal  to  the  dissipation  rate  at  viscous 
scales.  Although  his  model  adequately  represents  the 
initial  roll-off  of  the  velocity  spectrum  at  small  scales, 
measured  spectra  at  Kolmogorov  scales  roll  off  more 
rapidly  than  Pao’s  model  predicts.  Saffman  [1963]  de¬ 
rives  a  form  of  the  spectrum,  for  k  >  k,,  where  k,  = 
is  the  Kolmogorov  wavenunober,  which  rolls 
off  more  rapidly  than  Pao’s  model  predicts.  Panchev 
and  Kesich  [1969],  noting  that  the  assumptions  un¬ 
derlying  Pao’s  model  cannot  hold  as  the  peak  of  the 
shear  spectrum  is  approached,  propose  an  “interpola¬ 
tion  form”  which  joins  Pao’s  spectrum  at  low  wavenum¬ 
bers  to  Saifman’s  spectrum  at  high  wavenumbers.  In 
nondimensional  form,  the  three-dimensional  spectrum 
is 
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Figure  Cl.  Comparison  of  the  Pao,  Nasmyth,  and 
Panchev  and  Kesich  high-wavenumber  shear  spectra. 
The  Panchev  and  Kesich  spectrum  is  plotted  with  two 
values  of  the  Kolmogorov  constant,  a  =  1.6  and  2.0. 


where  a  is  the  Kolmogorov  constant  and  i/  is  the  kine¬ 
matic  viscosity.  To  compare  this  with  our  shear  probe 
measurements,  we  convert  to  the  one-dimensional  trans¬ 
verse  velocity  form  using 


EiCkx)  _  1 

(ci/S)!/'*  2 


which,  upon  a  change  of  variables  to  <  =  ki/k  where 
k  =  k/k,  and  fci  denotes  a  component  of  the  wavevector, 
yields 


ExCki) 

(f„5)l/4 


(4  (!)’)<  (C2) 

Plotted  in  Figure  Cl  are  the  Pao,  Nasmyth,  and 
Panchev  and  Kesich  spectra,  with  a  =  16  for  the  two 
theoretical  spectra.  Although  all  three  asymptotically 
approach  the  same  curve  at  low  wavenumbers,  their  be¬ 
havior  near  the  peak  of  the  shear  spectrum  and  in  the 
roll-off  differs  markedly.  The  Pao  spectrum  rolls  off 
most  gently,  whereas  the  Nasmyth  spectrum  is  the  most 
peaked.  The  shape  of  the  high  wavenumber  roll-off  is 
similar  for  the  two  theoretical  spectra;  the  Nasmyth 
spectrum  is  almost  linear  on  the  log-log  plot,  like  a 
power  law. 

Spectra  from  the  first  billows  sampled  are  better  fit  by 
the  Panchev  and  Kesich  spectrum  with  a  =  2.  The  ef¬ 
fect  of  changing  this  constant  is  illustrated  in  Figure  Cl; 


the  spectral  level  m  the  inertial  subrange  is  increased, 
and  the  wavenumber  where  the  spectrum  begins  to  roll 
off  is  decreased 
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